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We calculate the stationary current through two double quantum dots which are interacting via 
a common phonon environment. Numerical and analytical solutions of a master equation in the 
stationary limit show that the current can be increased as well as decreased due to a dissipation 
mediated interaction. This effect is closely related to collective, spontaneous emission of phonons 
(Dicke super- and subradiance effect), and the generation of a 'cross-coherence' with entanglement 
of charges in singlet or triplet states between the dots. Furthermore, we discuss an inelastic 'current 
switch' mechanism by which one double dot controls the current of the other. 
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I. INTRODUCTION 

The interaction with a dissipative environment can 
considerably modify the physics of very small systems 
which are described by a few quantum mechanical states 
only. 1 One may think of an excited atom that decays via 
the emission of a photon due to the coupling to the ra- 
diation field^ The influence of the environment becomes 
even more important if not only a single system is cou- 
pled to it but many. This introduces an indirect interac- 
tion between the otherwise independent systems which 
can result in an entanglement and collective effects of 
the small systems. In the case of identical excited atoms, 
the interaction to the common radiation field strongly 
affects the emission characteristics and leads to a collec- 
tive spontaneous emission, the so-called superradiance, 
as first pointed out by R. H. Dick o 3 i 4 i 5 nearly half a cen- 
tury ago. 

The influence of a dissipative environment on a single 
two-level system, the smallest non-trivial quantum sys- 
tem, has been studied extensively with the spin-boson- 
model- where the environment is modeled by a contin- 
uum of harmonic oscillators. Especially useful for the 
experimental realization of two level systems are cou- 
pled semiconductor quantum dots as these allow tun- 
ing of the parameters over a wide range£L2i2ii& More- 
over, in these systems transport spectroscopy is possible 
by connection with leado 8 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 : 22 . 
The dissipative environment is given by the phonons of 
the sample and governs the inelastic current through 
the system, 8 . 10 : 23 ! 24 ! 25 ! 26 . 27 : 28 The electron spin 29 ! 30 ! 31 : 32 
or the electron charg o 22 . 33 in quantum dots have also 
been suggested to provide a controllable realization of 
scalable qubitsi 29 i 30 i 34 i 35 i 36 i 37 Arrays of double quantum 
dotsS correspond to charge qubit 'registers', and sim- 
ple 'toy' models of N coupled two-level systems have 
been used to study collective decoherence effects in qubit 
registers Furthermore, controllable two-level sys- 
tems with Cooper pairs tunneling to and from a super- 
conducting box have been realized experimentally^S*^ 3 - 

Coherent effects in small clusters of two level systems 



caused by the coupling to a common environment have 
been realized mainly in the field of quantum optics. In 
ion laser traps, Dicke sub- and superradiance has been 
measured by DeVoe and Brewer— in the spontaneous 
emission rate of photons from two ions as a function of the 
ion-ion distance. Furthermore, entanglement in linear 
ion-traps can be generated by the coupling of (few-level) 
ions to a common single bosonic mode, the center-of-mass 
oscillator (vibration) mode 4 ^ 6 ^ 7 -. Even the generation 
of entangled light from white noised has been suggested. 

The appearance of collective quantum optical effects in 
mesoscopic transport has re-gained considerable interest 
quite recently. Shahbazyan and RaikbJ 9 - first predicted 
the Dicke (spectral function) effect to appear in reso- 
nant tunneling through two impurities, which was later 
generalized to scattering properties in a strong magnetic 
field 5 ^. The Dicke effect was predicted theoretically in 
'pumped', transient superradiance of quantum dot arrays 
coupled to electron reservoirs^, and in the AC conduc- 
tivity of dirty multi-channel quantum wires in a strong 
magnetic field i 5 ^ 

In this work, we focus on coherent effects in meso- 
scopic few level systems. As a realization, we choose 
two nearby but otherwise independent double quantum 
dots coupled to the same phonon environment. We study 
the influence of the resulting indirect interaction on the 
transport properties and calculate the stationary cur- 
rent. Signatures of 'super'- and sub-radiance of phonons 
are predicted which show up as an increase or a de- 
crease of the stationary electron current. We demon- 
strate that this effect is directly related to the creation of 
charge wave function entanglement between the two dou- 
ble dots, which appears in a preferred formation of either 
a (charge) triplet or singlet configuration, depending on 
the internal level splittings and/or the tunnel couplings 
to the external electron leads in both sub-systems. Gen- 
eration of entanglement via phonons becomes attractive 
in the light of recent investigations of single-electron tun- 
neling through individual molecules 5 ^ 5 ^ 6 * 5 ^ 5 ^, or quan- 
tum dots in freestanding 5 ?i 60 i 61 i 62 and movablo 63 i 64 i 65 i 66 
nano-structures, in both of which vibration properties on 
the nanoscale seem to play a big role. 
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The outline of the paper is as follows. In section[H]we Introducing the operators 
introduce the model and our method. Current superradi- 
ance is discussed in section ITTT1 Section Hvl presents cur- 
rent subradiance and the inelastic current 'switch' mech 
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anism. Finally, we conclude in section F7\ p t = \L,i) (R,i\ p\ = \R,i) (L,i\ (1) 

sl,i = |0, i) (L, i\ s Ryi = |0, i) (R, i\ , 



II. MODEL AND METHOD 

Our model is a system ('register') of two double quan- 
tum dots (DQDs), each of which consist of two individual 
quantum dots (called 'left' and 'right' in the following). 
Both double dots are coupled to independent left and 
right leads as depicted in Fig. \T\A. 

We concentrate on boson-mediated collective effects 
between the DQDs originating from the coupling of the 
whole system to a common dissipative, bosonic bath that 
will be specified below. In the following we completely 
neglect static tunnel coupling between the individual 
DQDs and, more important, inter-DQD Coulomb corre- 
lations. Although this is a severe limitation for the gen- 
eral applicability of the model, it still grasps the essential 
physics of dissipation induced entanglement. However, 
one might envisage configurations with intradot Coulomb 
matrix elements much larger than interdot matrix ele- 
ments. 

In this paper, we choose the simplest possible de- 
scription of an environment coupling in close analogy 
to the standard spin-boson Hamiltoniani. The results 
of this model for the tunnel current through one double 
dot are in relatively good agreement with experimental 
observations 2 ^ 3 -. The role of off-diagonal terms in a sin- 
gle DQD has been discussed recently 67 . 



Hamiltonian 



the total Hamiltonian can be written as 

N r 

H = ^ e L ,i n L ,i + en,i n R:i + T C:i (pi + p\) 

i 

+ XI y M c li s L,i + h.c. + 4,i c M 
k fc 

+ V * d b s R,* + h - c - + I>f 4, A* 
i i 

+ X! 7 <?- 1 ( a <3 + a -q)( n L,i ~ n R.i) +^Lu q a\a q 



(2) 



Here, the electrons in mode fc (/) with energy e^mj m 
the left (right) leads pertaining to DQD i are described 
by creation operators cl i (dj f ), and the coupling matrix 

elements to the leads are denoted by V k / .A boson in 
mode q with energy uj q is created by the operator a q . As 
in the standard spin-boson model, we assume a simplified 
coupling to the quantum dots which is purely diagonal 
with matrix clement 7 9 j for mode q to the ?-th double 
dot. 

So far, no further assumptions have been made with 
respect to the specific realization of the DQDs and the 
dissipative bath. Nevertheless, the system we have in 
mind are lateral or vertical double dots, where the pri- 
mary bosonic coupling has been shown due to phonons 
of the semiconductor substrate. The microscopic details 



determine the tunnel matrix elements T r , , V, 



L/R 



the electron-phonon coupling constants 7 g) j. 



k . i 



and 



The Hamiltonian and the subsequent derivation of the 
master equation is given for the general case of N double 
quantum dots. We study the stationary tunnel current 
through the dots with all lead chemical potentials such 
that electrons can only flow from the left to the right. 
Furthermore, we restrict ourselves to the strong Coulomb 
blockade regime in each individual double dot where only 
one additional electron is allowed on either the left or 
the right dot. The Hilbert space of the i-th double dot 
then is spanned by the three many-body states \L, i) (one 
additional electron in the i-th left dot at energy £l,i), 
\R,i) (one additional electron in the z-th right dot at 
energy £r,j), and |0,i) (no additional electron in cither 
of the dots) . For N = 1 this has been proveniSiSi to be a 
valid description of non-linear transport experiments in 
double quantum dots&iS. 



B. Density matrix 

In the following, we employ a master equation de- 
scription for the time evolution of the register within 
the Born-Markov approximation, which takes into ac- 
count the interactions with the leads and the bosonic 
environment up to second order. Alternatively, electron- 
phonon interactions can be treated exactly by a polaron 
transformation 22 23 and perturbatively in the tunnel ma- 
trix elements T c ^. For T C) j < \en — £Ri\ and small cou- 
pling to the bosonic bath, the results of both methods 
practically coincide^. 

The time derivative of the reduced density matrix p(t) 
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FIG. 1: A: N — 2 'charge qubit register' with two double 
quantum dots coupled to independent electron leads. B: en- 
ergy diagram of one individual double dot. 



of the double quantum dots is given by 

p(t) = - J dt' (Tr Ros , c { [H c (t), [H e (t% p{t) ® Ro >e ]] } 

+Tr Rcs , p { [H p (t), [H p (t'),p{t) ® flo, p ]] } ) , 

(3) 

where the tilde indicates the interaction picture, H e (H p ) 
denotes the interactions between the double dots and the 
leads (the phonons), and Rq. c {Ro,p) is the density matrix 
of the leads (the phonons) . Equation J2J is the sum of an 
electron and a phonon part since we neglect correlations 
between leads and phonons. 

The trace over the equilibrium electron reservoirs, 
TrRes.e, results in Fermi functions of the leads. As we 
are interested in large source-drain voltages between the 
left and the right leads, the Fermi functions of the left 
leads can be set to one and those of the right leads to zero. 
Moreover, the energy dependence of the tunnel rates 



L/R,i 



_L/R, 



(4) 



is neglected. 



and the geometry of the phonon substrate (bulk, slab 27 , 
sheet etc.), the ^ q ^i will never be exactly identical in real 
situations. Therefore, Eq. JSJ can only be regarded as an 
idealized limit of, e.g., a phonon resonator or a situation 
where the distance between different double dots is small 
as compared to the relevant phonon wavelengths. 

We define a correlation function of the boson system 



K(t) 



du> p(to] 



1 - e-P" 



(6) 



that results from the trace over the bosonic degrees of 
freedom. Here, (3 = 1/ksT denotes the inverse phonon 
bath temperature, and the spectral function p(oj) of the 
bosonic environment is defined as 



P{w) = ^l7 9 | 2 5{u~uj q ). 



(7) 



For the calculations, we use the spectral function of bulk 
acoustic phonons with piezoelectric interaction to elec- 
trons in lateral quantum dots 2 ^ 2 ^, 



p(uj) = gu) ( 1 sin ( — | ) < 



(8) 



where g is the dimensionless interaction strength, u> c the 
cut-off frequency and the frequency u)d is determined by 
the ratio of the the sound velocity to the distance between 
two quantum dots. 

In the following, integrals over K{t) are required as 

T c ,i = J Q K(t) cos(A^) dt = | p(Ai) coth(^), 

T s ,i = / K{t) sin(A l t) dt =-i- p(Ai). 
Jo z 

(9) 



with the hybridization energy = (ef + 4T ( 2 i ) 1 / 2 and 
the energy bias £i = £l,i — £r.i in the z-th dot. The in- 
tegrals are calculated neglecting the principal values 2 ^. 
We furthermore assume a spectral function p(A 4 ) such 
that Tc,i — * for A,; — > which is fulfilled for micro- 
scopic models of the electron-phonon interaction in dou- 
ble quantum dots 221 " ' 



Electron-phonon interaction 



In the following, we consider identical electron-phonon 
interaction in the DQDs, 

Iq-i = lq ■ (5) 

Depending on the relative position of the quantum dots 
(lateral, vertical), the electron wave functions in the dots, 



2. Master equation 

Inserting the traces over the electron reservoirs and 
the bosonic bath into Eq. © and transforming back to 
Schrodinger picture yields a master equation for the re- 
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duced density matrix of the total DQD register, 

N 



p(t) = i ^2 | \p(t) , £l, 4 n L .i + £ Rii n R .i + T c .i(p t + p\ 

i=l 

+ -y 1 ( 2 4,i />(*) S i,i - S £,* S L ~ S ^.i S L,i) 

+ ( 2 P{t) s R , % ~ s R<i s R . t p(t) - p(t) s R4 s R .^ I 
- ^ { \jnL,i-n Rt j), Aj p(t) - (n L ,i-n Rji ),p(t)A^ |, 

(10) 

with 



(11) 



From Eq. @ it is obvious that the influence of the 
bosonic bath enters only via the spectral functions p(ui) 
as defined in Eq. (J7J- All microscopic properties of the 
phonons and their interaction mechanism to the electrons 
in the quantum dots are described by these functions. 

Furthermore, we point out that the mixed terms i =/= j 
in Eq. (|10|l are responsible for the collective effects to 
be discussed in the following. Without these terms, the 
master equation would merely describe an ensemble of N 
independent DQDs. In that case, an initially factorized 
density matrix of the total system would always remain 
factorized and no correlations could build up. The terms 
i =/= j introduce correlations between the different double 
dots, the origin of which lies in the coupling to the same 
bosonic environment. 



III. CURRENT SUPERRADIANCE 

We restrict ourselves to the stationary case where the 
time derivative of the density matrix, p(t), vanishes. 
Then, Eq. I|10(l reduces to a linear system of equations 
which can be easily solved numerically. Results for a 
single double quantum dot, N = 1, can be obtained 
analytically^ and are given for two expectation values 
below, Eq. (|T7jl . For N > 1, the dimension of the density 
matrix grows as 9 N (although not all of the matrix ele- 
ments are required) whence analytical solutions become 
very cumbersome. For the rest of this paper, we restrict 
ourselves to the case of two double dots (N = 2), called 
DQD 1 and DQD 2 in the following. 
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FIG. 2: Total current through two double quantum dots as a 
function of the bias £1. The parameters are T c ,i ~T C: 2 = 3/ieV, 
r.6,1 =Fr > i = Fl,2 =rH,2 = 0.15/ieV, and for the spectral func- 
tion g = 0.01, T = 23mK, u d = 10/xeV and lo c = lmeV. These 
values are used throughout the whole article if not stated oth- 
erwise. 



cannot tunnel between different double dots. The current 
operator of DQD i is 



(12) 



and the corresponding expectation values can be ex- 
pressed by the elements of the density matrix as 



h — r— ImjpiiiLL + Prrlr + Porlo}, 



h = - 



h 

2 T c , 2 e 



I m {PRLLL + PRRRL + PROOl}, 



(13) 



with the notation 



Pjii'j'= 2 01 ® P ® \j') 2 > i,je{L,R,0}. 

(14) 

The set of linear equations corresponding to Eq. IjlUfl for 
N = 2 is given in appendix^ Eq. (|A1|) . From the numer- 
ical solution of Eq. IjAlfl we find the stationary current 
through two double quantum dots as a function of the 
bias £1 in the first double dot while the bias £2 in the 
second is kept constant, cf. Fig. [21 The overall shape 
of the current is very similar to the case of one individ- 
ual double quantum doljSii 2 ^, with its strong elastic peak 
around £1 = and a broad inelastic shoulder for £1 > 0. 
The interesting new feature here is the peak at the res- 
onance £1 = £2 which is due to collective effects to be 
analyzed now. 



A. Stationary current 



Cross coherences 



The total electron current is simply given by the sum 
of the currents through the individual DQDs, as electrons 



The effective interaction between the two DQDs re- 
sults from the simultaneous coupling of both double dots 
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to the same phonon environment. It appears in the mas- 
ter equation l|l(J|l as the mixed terms i ^ j in the sum. 
In the explicit form of the master equation IjAljl . the ef- 
fective interaction is connected to six matrix elements 
only (and their complex conjugates). These elements are 
Prlll, Plrll, Prrlr and prrrl, all of which enter the 
expression for the current, Eq. (| 1 3|) . and the two 'cross 
coherence' matrix elements 

PRLRL = (PiPs) , PRRLL = (PIP2)- (15) 

Therefore, we approximate the collective effects caused 
by the effective interaction starting from the solution of 
the non-interacting master equation, without the mixed 
terms i^j, and assume that only those matrix elements 
mentioned above are affected by the interaction. 

In the non-interacting case, the cross coherence is sim- 
ply the product of the corresponding matrix elements of 
independent double dots, 

(pIp*) = (p\) (P2) , (pm) = (pi) (p 2 ) . (16) 

These can be solved analytically, 

+ T RJ (iT cd + ij){ie 3 + -T RJ + 2a,)) , 

(n L ,j) = 1 - 7: ' jT, ; w '' " /j (2 ,,, + T^Tnj + 4a,)) , 

(17) 

where (tilj) is given for later reference, ay, 0j, and "fj 
as defined in the appendix, Eq. (|A3|> . and with 

+ r c 2 J (2r Lj +r^)(r^+4a,) + 2t c ^ j1j (t l , +t r ^. 

(18) 

In the inelastic regime, T c <C e, of the non- interacting 
case, the cross coherences (p\p2) and (P1P2) tend to zero 
as can be seen from Eq. H17|) . Moreover, we neglect the 
imaginary part of the cross coherences in the interacting 
case. Then, the change in the current through DQD 1 
due to collective effects can be approximated by 

Ah = (Re{(p\p 2 )}-Rc{( PlP2 )}). (19) 

Correspondingly, the change AI 2 of the current through 
the second double dot DQD 2 is obtained from AI± by 
exchanging the subscripts 1 and 2. Hence, the alteration 
in the current is proportional to the real parts of the cross 
coherences (p\p2) and (pip 2 ) between the two DQDs, 
which confirms the collective character of the effect. This 
result is corroborated by plotting the real parts of the 
cross coherences as a function of e 1; cf. Fig. [3] One 
recognizes that (p\p 2 ) is peaked around £1=62, whereas 



{P1P2) has a peak at E\ = —£2- The increase of the current 
at £1 = £2 is therefore due to the maximum of the first 
correlation (p\p2)- 

If we neglect the changes of all other elements of the 
density matrix that are caused by the effective interac- 
tion between the two DQDs, the real part of the cross 
coherence (p\p 2 ) can be approximated around the reso- 
nance £1 =£2 as 




\ (r_R,i + r fij2 ) 
(Ei^F + KiV + r^)^ 



^7iRe{(p 2 )} (nL,i) + 72Re{(pi)} (n^.a))- 

(20) 

One recognizes that (jp\p2) is Lorentzian shaped as a 
function of the energy difference £1— £2- The result of 
Eq. (121 III with (pj) and (riL.j) as given in Eq. (|17|1 is in 
good agreement with the numerical solution of the master 
equation lfTU|) (inset of Fig. |2J. 

Next, we insert the result for the cross coherence in 
Eq. I|19|) and find for the change of the tunnel current due 
to interaction effects between the two double quantum 
dots around the resonance £1 =£2: 

_ eT c .i 72 Pr,i+Tr,2 

1_ hex (Ei-^ + ^r^+r^^ (21) 

(71 Re{ (p 2 ) } (ni,i) + 72 Re{ (pi) } {n L , 2 ) ) • 

Again, the change in the current through the second dou- 
ble dot, A/2, is obtained by exchanging the subscripts. 
This approximation overestimates the actual change in 
the current for the parameters chosen in the previous 
section but provides a good qualitative description for 
the effect of the enhanced tunnel current. A comparison 
between this result and the numerical solution is given 
below. 



C. Singlet and triplet states 

The collective effects in the two double quantum dots 
are connected with the cross coherence function (p\p2), 
Eq. (|15fl . For the 'two-qubit register' one can easily prove 
the operator identity 

pIp2+pIpi = Pt -P So , (22) 

where P^, is the projection operator on the state \ip), 
P.0 = (ip\, and triplet and singlet do not refer to the 
real electron spin but to the 'pseudo' spin defined in the 
two dimensional Hilbert space span(|L), \R)), 

|T+) = |i> 1 |L) 2l |T_) = |i?) 1 |i?) 2 , 

\T ) = ^=\[\L) 1 \R) 2 + \R) 1 \L) 2 ), (23) 

|5o) = ^(|i) 1 |-R) 2 -|ii) 1 |i) 2 ). 



6 



0.08 



0.06 - 



0.04 



0.02 




FIG. 3: Real parts of the cross coherences from the master 
equation ffllOl) as functions of the bias in the first double dot 
(s2 =30/xeV and the other parameters agree with Fig.|2J. The 
inset compares the approximation for Re{(pjp2)}, Eq. I20H 
(dotted line), with the solution of the master eq. (solid line). 



With 2Re( P \p 2 ) = (Pro) - (Psr ) and the proportional- 
ity All oc Re(p{p 2 ) for E\ w e 2 , cf. Eq. (jT§)) . it follows 
that the current enhancement A/i is due to an increased 
probability of finding the two electrons in a (pseudo) 
triplet rather than in a (pseudo) singlet state. In the 
following, we demonstrate that the mechanism underly- 
ing this effect is indeed the Dicke superradiance effect 
known from quantum optics. 



D. Dicke effect 

Superradiance emerges in the collective spontaneous 
emission from an ensemble of identical two-level atoms. 
If N excited atoms are concentrated in a region smaller 
than the wavelength of the emitted radiation, they do 
not decay independently anymore. Instead, the radiation 
has a higher intensity and takes place in a shorter time 
interval than for an ensemble of independent atoms due 
to the coupling of all atoms to the common radiation 
field. 

Let us now consider the case N — 2 and calculate (simi- 
lar to the original work of Dicke^) the decay rate T of two 
initially excited atoms with dipole moments d\ and d 2 ) 
at position r x and r 2 due to the interaction with light, 



H ep h — ^2 9Q ( a -Q 



Q 



(24) 



x \di expi(Qr 1 ) + d 2 expi(Qr 2 ) 
from which the spontaneous emission rate of photons 



with wave vector Q follows (Fermi's Golden Rule), 

r±(Q) \9q\ 2 \ 1 ± ex P [*Q( r 2 - n)]l 2 <^o - uq), 

Q 

(25) 

where Q = ujq/c, ujq is the transition frequency between 
the upper and lower level, and c denotes the speed of 
light. The interference of the two interaction contri- 
butions die 1 ^^ 1 ' and e^e 1 ^ 1 " 2 -' leads to a splitting of 
the spontaneous decay into a fast, 'superradiant', decay 
channel (T + (Q)), and a slow, 'subradiant' decay channel 
(r_(Q)). This splitting is called 'Dicke-effect'. 

Loosely speaking, the two signs ± correspond to the 
two different relative orientations of the dipole moments 
of the two atoms. More precisely, from the four possi- 
ble states in the Hilbert space of two two-level systems, 
H9 = C 2 <X> C 2 , one can form singlet and triplet states 
according to \S ) := ^ (| U> - UT», ' T +> := ' TT >' 
\T ) := ^ (| U) + I IT)), and |T_) := | U). The super- 
radiant decay channel occurs via the triplet and the sub- 
radiant decay via the singlet state o?'^? . In the extreme 
'Dicke' limit where the second phase factor is close to 
unity, exp [iQ(r 2 — ri)] 1, it follows that r^(Q) = 
and T + (Q) = 2T(Q) where T(Q) is the decay rate of 
one single atom. This limit is theoretically achieved if 
|Q(r 2 — rx)| <C 1 for all wave vectors Q, i.e. the distance 
between the two atoms is much smaller than the wave 
length of the light. 

We mention that in practice, this 'pure' limit, where 
the subradiant rate is zero and the superradiant rate 
is just twice the rate for an individual atom, is never 
reached. In a recent experimental realization of sub- and 
superradiance from two laser-trapped ions, DeVoe and 
Brewer— measured the spontaneous emission rate of pho- 
tons as a function of the ion- ion distance in a laser trap of 
planar geometry which was strong enough to bring the 
ions (Bajjg) to a distance of the order of 1/um of each 
other. 

The two double quantum dots behave in analogy to 
the two atoms considered above. For a positive bias, 
e = el — £r > 0, the state \L) can be identified with the 
excited state and \R) with the ground state. The inelastic 
rate v with which \ L) decays to \R) can be calculated with 
Fermi's Golden Rule, 



8-kT 2 



-4T2) 



(26) 



In contrast to a two level atom, a third state |0) exists 
in the double quantum dot as the additional electron can 
tunnel into the leads. 

The use of triplet and singlet states as defined in 
Eq. 123fl allows us to find an analytical result for the 
stationary current that quantitatively coincides with the 
exact numerical solution extremely well. We consider the 
rate equation for the probabilities of the corresponding 
nine states and take into account the doubling of the in- 
elastic rates due to the Dicke effect in the triplet channel, 
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FIG. 4: Enhancement of the tunnel current A/i at the res- 
onance £i = Si = 30/ieV as a function of the dimensionless 
electron phonon coupling constant g, Eq. ©. The additional 
current vanishes at <7~0.02 when the tunnel rates to the dou- 
ble dot and between the dots become equal, v = Y. The inset 
shows the difference in probabilities for triplet and singlet. 



The difference All = 1\ — l\ represents the additional 
current due to the Dicke effect and is shown in Fig.0]as a 
function of the dimensionless coupling strength g to the 
bosonic environment, together with a comparison to the 
All as obtained from the numerical solution of Eq. l|10fl . 
Both results agree very well, indicating that it is indeed 
the Dicke effect that leads to the increase in the tunnel 
current. In addition, we show (inset of Fig. 0} the differ- 
ence between triplet and singlet occupation probability 
that follow from the Eq. (|27Jl as 



PT - PSo 



2x(x + 2)(x- 1) 
9a; 3 + 23a; 2 -t- llx + 2' 



(30) 



This is in excellent agreement with the numerical results 
and underlines that the change in the tunnel current due 
to collective effects is proportional to px —Ps > as already 
discussed above. This demonstrates that the effect of 
superradiance amplifies the tunneling of electrons from 
the left to the right dots resulting in an enhanced current 
through the two double quantum dots. 



IV. CURRENT SUBRADIANCE AND 
INELASTIC SWITCH 



Poo 
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Pol 

Pro 
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PT+ 
PT 

i>T~ 

PSo 



1 



T PT 
T PT 



2 

Tp T 



2I>oo, 
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+ 2 
1 
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^Pro, 



fps 
FPSO 



Tpoo - 

Tpoo - 

vPlo 
vpoh 

I>0L 

2vp T + 
2vp To - 

2 r PoR + 2 r PR0 - Fpso ■ 



{v 



Tp T _ - 2Tpo R , 
-Tp L0 - 2vp T+ , 



(27) 



1 



1 



- ^poR + - rp m - (2^ 

2Fp T _, 



Here, identical tunnel rates to all four leads have been 
assumed, ^l,x = ^r,x=^l,2 = ^r > 2 = ^', and plq denotes 
the probability to find the first double dot in state \ L) and 
the second in state |0). Electrons can also tunnel into and 
out of the singlet state due to the coupling to the leads 
which is not possible in the original Dicke model. In the 
stationary case, the Eq. I|27fl can be easily solved. For the 
current through one of the two double dots we obtain 



h 



eT 

~h 
eT 

IT 



(Poo + Pol 

x(Ax + 1) 
9x 2 + 5x + 1 



Por) 

x = v/T. 



(28) 



This can be compared with the tunnel current through 
one independent double dot, 7° obtained by a similar rate 
equation, 

' ' ' (29) 



J? 



h l + 2x 



The close analogy with the Dicke effect suggests the 
existence of not only current super-, but also current 
subradiance in the register. In the subradiant regime, 
the two DQDs form a singlet state where the tunneling 
from the left to the right quantum dots is diminished, 
resulting in a weaker tunnel current through the dots. 



A. Current antiresonance 

Subradiance occurs in our system in a slightly changed 
set-up where electrons in the second double dot are pre- 
vented from tunneling into the right lead, = 0, as 
indicated in the inset of Fig. |3J Then, the additional 
electron is trapped and no current can flow through the 
second double dot. Nevertheless, this electron can affect 
the tunnel current through the first double dot: Instead 
of a maximum, we now find a minimum at the resonance 
£\ =£2- Fig.[S]shows how the positive peak in the current 
I\ develops into a minimum as the tunneling rate Tr^ is 
decreased to zero. This minimum is indeed related to an 
increased probability of finding the two dots in the sin- 
glet state \Sq) rather than in the triplet state \Tq), as can 
be seen from the inset of Fig.JSJ Thus, in this regime the 
effect of subradiance dominates, leading to a decreased 
current. 

This behavior is again consistent with the approxima- 
tion Eq. (|20|) for the cross coherence (jp\p2)- Taking into 
account the different non-interacting matrix elements in 
the two double dots, ( n L,2) and (pi) ^ (P2) due 

to Tr,i 7^ Yr.2, we find a negative cross coherence at 
the resonance from Eq. This corresponds to an in- 

creased probability for the singlet state and according 
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ei/fieV Ei/fieV 



FIG. 5: Transition from an increased to a decreased current 
through the first double quantum dot for different tunnel rates 
Ti?,2 (in neV) and £2 = 30^eV. The left inset shows schemat- 
ically the set-up for I\r,2 = and the right inset gives the 
difference of triplet and singlet for the same case. 



to Eq. 121fl to a negative peak in the tunnel current, in 
agreement with our numerical solution. 



B. Inelastic current switch 

Up to now, we have regarded the cross coherence 
{p\p2) and its effects on the current only at the reso- 
nance £i = £2- However, it was already pointed out in 
section IlII Bl that another cross coherence, (pi^), ex- 
hibits a resonance if the bias in one dot equals the neg- 
ative bias in the other dot, E\ — —Ei (cp. Fig. |3J). This 
case is considered in the following. 

We use a fixed negative bias £2 < in the second dou- 
ble dot as indicated in the inset of Fig. [S] Consequently, 
electrons cannot tunnel from the left to the right dot such 
that the second double dot is blocked and no current can 
flow through it. The presence of the first double dot, 
though, lifts this blockade and enables a current through 
the second double dot if the resonance condition E\ = —£2 
is fulfilled. The current I2 is shown in Fig. as a func- 
tion of the bias in the first double dot, E\. Due to the 
coupling to the common phonon environment, energy is 
transferred from the first to the second double dot, allow- 
ing electrons to tunnel from the left to the right in the 
second double dot. At the same time, the current through 
the first double dot is decreased (not shown here). 

We can approximate the current through the second 
double dot around £1 = —£2 taking into account only 
(P1P2) in Eq. (|19fl . A similar calculation as for (jp\p2), 



FIG. 6: Tunnel current through the second double dot which 
is blocked due to a negative bias, £2 = — 50/xeV, as depicted 
in the inset (g = 0.015). The approximation for A/2, Eq. 
(dashed line) , agrees well with the result of the master equa- 
tion 1101 (solid line - the finite offset of which is the tail of 
the elastic current at £2=0). 

Eq. (201), gives 

_ 2T C .2 71 e ^R,i + ^R,2 + 8a 

2 ~ £2h (£i+£ 2 ) 2 + (ir flil +ir fli2 +8a) 2 

( 7l Rc{(p 2 )} (n L ,i) +7 2 Re{(pi)} (n L , 2 )), 

(31) 

with a — a\~ci2 evaluated at the resonance, where both 
systems are identical except of the bias. This approxi- 
mation again is in good agreement with the numerical 
solution of Eq. (ID) , as can be seen from Fig. ED 

Our results suggest that the current through one of the 
DQDs can be switched on and off by appropriate manip- 
ulation of the other one. We emphasize that this mecha- 
nism is mediated by the dissipative phonon environment 
and not the Coulomb interaction between the charges. 
As this effect is very sensitive to the energy bias, it al- 
lows to detect a certain energy bias in one double dot by 
observing the current through the other double dot. 

V. CONCLUSION 

In this work, we have investigated collective effects in 
two double quantum dots. An indirect interaction arises 
between the two double dots due to the coupling to the 
same phonon environment. We predict that the Dickc 
effect causes a considerable increase or decrease of the 
tunnel current, depending on the choice of the parame- 
ters. The occurrence of the Dickc effect in the transport 
through mesoscopic systems has already been pointed out 
by Shahbazyan and Raikk 4 ^. In their system, the cou- 
pling to the same lead is responsible for collective effects. 
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Usually, the Dicke effect manifests itself in a dynamic 
process like the spontaneous emission of an ensemble of 
identical atoms 44 - 69 . Transport through double quantum 
dots, however, allows to study a time independent form 
of the Dicke effect. Moreover, we have demonstrated that 
the change of the tunnel current is connected with an en- 
tanglement of the different double dots. This opens the 
possibility to realize and to measure specific entangled 
states of two double dots. In particular, one can switch 
from a predominant triplet superposition of the two dou- 
ble dots connected with an increased tunnel current to a 
predominate singlet state leading to a reduced current. 

The results discussed here were derived for the ideal 
case of an identical electron-phonon coupling in both 
double quantum dots. Furthermore, the Coulomb inter- 
action between the two double dots has not been consid- 
ered here. In a real experiment, these assumption will 
never be perfectly fulfilled and would lead to deviations 
from the collective effects presented above. However, 
we predict that even in presence of inter-dot Coulomb 
interactions, phonon mediated collective effects should 
persist as long as a description of the register in terms 
of few many-body states is possible. These many-body 
states (that would depend on the specific geometry of 
the register) would than replace the many-body basis 
{|0, i), \L, i), \R, i)} (i — 1, 2) used in our model here. 

We have derived the master equation for the general 
case of N double dots but only focused on N = 2 which 
is the simplest case where collective effects occur. In gen- 
eral, one of the main characteristic features of superradi- 
ance is the quadratic increase of the effect with increasing 
number of coupled systems. For the spontaneous collec- 
tive emission from N excited two level atoms, this means 
that the maximum of the intensity of the emitted radi- 
ation increases with the square number of systems, ./V 2 , 
while the time in which the decay takes place decreases 
inversely to the number of systems, 1/N. Therefore, we 
expect that the collective effects as presented here be- 



come even more pronounced if more than two double dots 
are indirectly coupled by the common phonons. 
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APPENDIX A: MASTER EQUATION FOR TWO 
DOUBLE QUANTUM DOTS 

The dimension of the density matrix p for N double 
quantum dots is equal to 9 N such that the master equa- 
tion IjlOjl corresponds to 81 coupled differential equations 
for N = 2. It is, however, not necessary to solve all 81 
equations as we study the current which requires the 
knowledge of only six matrix elements, cp. Eq. JTj5J. 
The smallest closed subset of equations, containing the 
equations for those six elements consists of 25 equations. 

The mixed terms in the master equation IjlOjl . i ^ j, 
describing the indirect interaction between the two DQDs 
due to the coupling to the same phonons, are marked in 
the following with an additional prefactor q. Setting q = 
results in the master equation for two completely inde- 
pendent double dots coupled to independent phonons. 
The interacting case corresponds to q = 1- Note that 
the elements of the density matrix are expressed with 
respect to the basis {\L) , \R) , |0)} for each double dot. 
Due to the tunneling of electrons between the left and 
right quantum dot, these states are no eigenstates of the 
unperturbed Hamiltonian. Finally, the master equation 
for the elements of the density matrix reads 



PLLLL — iTc,l{PLLRL " PLRLl) + iT c 2 {pLLLR — PRLLl) + ^L,l pLOOL + r_L,2 POLLO, 

PLLLR. = iT c ,l(pLLRR — PLRLr) + iT c ,l{PLLLL — PRLLr) + ^L,l PLOOR — 72 PLLLL — Pi PRLLR 

- (i£2 + 2^ R ' 2 + 2a2 ) PLLLR + q Pi (PLLRR - PLRLr) , 

PLLRL — iTc,\{PLLLL — PLRRl) + iT c , 2 {pLLRR — PRLRh) + ^L,2 POLRO — 7l PLLLL — Pi PLRRL 

~ (i£l + ^ r -R,l + 2a l) PLLRL + Q Pi (PLLRR - PRLRl), 
PLLRR — iTc,\{PLLLR — PLRRr) + iT c2 (fiLLRL — PRLRr) — 7l PLLLR — Pi PLRRR 

- 72 PLLRL ~ Pi PRLRR - (l£l + %£l + ^ ' R.l + R.2 + 2oi + 2a 2 ) PLLRR 

q (2(ai + a2) PLLRR + Pi PRLRR + Pi PLRRR + 72 PLLRL + 7l PLLLr) , 
PRLLR — iT c ,l(pRLRR ~ PRRLr) + iT c ,l(PRLLL — PLLLr) + ^L,l PROOR — F R.2 PRLLR, 
PRLRL — iT c ,l(pRLLL — PRRRL) + iT c ,l{pRLRR — Pllrl) — 71 Prlll — Pi Prrrl 

1 1 

- 72 PLLRL ~ Pi PRLRR - (i£l~ %£l + 1^R,1 + ^ 'r,1 + 2ai + 2a 2 ) PRLRL 
+ q (2(oi +a 2 ) PRLRL + 72 PLLRL + Pi PRRRL + P2 PRLRR + 7l PRLLl) , 

Prlrr — iT c ,i{PRLLR — Prrrr) + iTc,2^.PRLRL — Pllrr) " 7i Prllr — Pi Prrrr 

- {i£\ + ~j^R,l + 2 + 2ai ) PRLRR + qi2{pLLRR ~ PRLRl), 
POLLO — iT Ct i(poLRO — P0RL0) + ^L,l P0000 + ^R,2 PRLLR — ^L,2 POLLO, 

POLRO — iTc,l{PQLLO — P0RR.0) + ^R,2 PRLRR. — 7l POLLO — Pi PORRO ~ + 1^ RA ^ L ' 2 ^ ai ) POLRO, 

PLRRL — iT Ct i(pLRLL — PLLRL) + iT c>2 (pLRRR ~ PRRRL) + ^L,2 PORRO — T R.l PLRRL, 

Plrrr — iT Ct i(pLjiLji — Pllrr) + iT c ,2(pLRRL — Prrrr) — 72 Plrrl — Pi Prrrr 

- (ie 2 + Tji^ + -j£R,2 + 202) PLRRR + qjl(pLLRR ~ PLRLr) , 

Prrrr — iTcs{prrlr — Prlrr) + iT c ^{pRRRL — Plrrr) — {^ra + ^R.,2) Prrrr, 

PORRO = iT c ,l{pvRL0 — POLRtt) + T ' R,2 PRRRR ~ R,l + ^L,l) PORRO, 
PLOOL = iT c ,2(pL00R ~ PROOl) + ^ R,l PLRRL + ^L,2 P0000 — ^L,l PLOOL, 

Ploor = iT Cj2 (pLooL — Proor) + T^a Plrrr — 72 Plool — Pi Proor — (i£2 + ^l,i + i^R,^ + ^02) Ploor, 

PROOR = iTc,2(pR00L ~ PLOOr) ~ ^L,l PROOR + ^R.l PRRRR ~ ^ R,2 PROOR, 
P0000 = ^R,l PORRO + Tfl : 2 PROOR — (^L,l + ^,2) P0000- 

I 

The remaining 8 equations follow immediately since p is 
an hermitian operator, 

Pjii'f = Pj>i>ij , (A2) 
and the coefficients aj, pj, and -fj are defined as 

4nT c i /A x , (P^i\ 
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